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In our first paper, we showed how a non-local effective Hamiltionian for short-ranged wetting 
may be derived from an underlying Landau-Ginzburg- Wilson model. Here, we combine the Green's 
function method with standard perturbation theory to determine the general diagrammatic form 
of the binding potential functional beyond the double-parabola approximation for the Landau- 
Ginzburg- Wilson bulk potential. The main influence of cubic and quartic interactions is simply to 
alter the coefficients of the double parabola-like zig-zag diagrams and also to introduce curvature and 
tube- interaction corrections (also represented diagrammatically) , which are of minor importance. 
Non-locality generates effective long-ranged many-body interfacial interactions due to the reflection 
of tube-like fluctuations from the wall. Alternative wall boundary conditions (with a surface field 
and enhancement) and the diagrammatic description of tricritical wetting are also discussed. 



I. INTRODUCTION 

In we showed how a non-local interfacial Hamil- 
tonian for short-ranged wetting 0] may be derived from 
a Landau-Ginzburg- Wilson model using a diagrammatic 
formalism based on Green's functions While the 

definition of the interfacial model is the same as that 
forwarded by other authors [1, its evaluation is 

non-perturbative in the interfacial gradient and reveals 
important non-local features. This has a number of ad- 
vantages 0, [sl over previous, local approximations, and 
appears to overcome a series ofproblems associated with 
short-ranged wetting I, i, i, i, S M M El, M M 
[Tsj . The interaction of the interface and wall is described 
by a binding potential functional W which has an elegant 
diagrammatic expansion 

W = ai^T^ +bi^\r + ... (1) 
— • ' — 

and an appealing physical interpretation as tube- like fluc- 
tuations [16j which zig-zag between the surfaces. 

In this paper, we demonstrate the robustness of the 
diagrammatic expansion by extending the derivation be- 
yond the double-parabola (DP) approximation for the 
bulk thermodynamic potential appearing in the Landau- 
Ginzburg- Wilson (LGW) Hamiltonian. We establish the 
general form of the asymptotic expansion of W, which 
now includes diagrams representing curvature corrections 
and tube-tube interactions. However, these are of negli- 



gible importance and the above diagrammatic expression 
remains valid albeit with numerical changes to the val- 
ues of the coefficients ai and &i, which can be determined 
exactly. 

Our article is arranged as follows: After recapping 
briefly the central results of [ij , sections IIIHIVII describe 
the detailed derivation of W beyond DP using the same 
diagrammatic formalism. The proof is rather technical 
and, to simplify things, we continue using flxed bound- 
ary conditions at the wall until the final section. Sections 
IVIimXI are a lengthy discussion of the interpretation of 
the model and the evaluation of the diagrams for wet- 
ting transitions at planar walls which extends the anal- 
ysis given in In particular, non- locality is shown to 
induce weak but long-ranged two-body interactions de- 
scribing the repulsion of the interface from the wall. Al- 
ternative forms of boundary conditions, including cou- 
pling to an external surface field and enhancement p^ . 
and also diagrams describing tricritical wetting are also 
discussed. 



II. WETTING DIAGRAMS WITHIN THE 
DOUBLE-PARABOLA APPROXIMATION 

We begin with some general considerations. Imagine a 
system bounded by a wall described by a height function 
'(/'(x) measured above some reference plane with parallel 
vector displacement x = (x, y) (see Fig.l). The wall is in 
contact with a fluid phase a but preferentially adsorbs a 
fluid phase /3 which forms a thin film intervening between 



A(|)(m) 




FIG. 1: Schematic diagram showing a wetting layer of phase P 

at a non-planar wall-a interface. Here, ^(x) and V(x) are the FIG. 2: The relative potential A(/>(m) in the "m''" theory Q 
interfacial collective coordinate and waU height, respectively. and DP approximation ^ at bulk coexistence. 



the bulk a phase and the wall. Throughout our paper, we 
will restrict our attention to the case of bulk two-phase 
coexistence, although it is easy to extend the analysis 
to non-zero values of the bulk ordering field. For sys- 
tems with short-ranged forces, a convenient microscopic 
starting pointing for studying this is the continuum LGW 
Hamiltonian 



H[m] 



dv 



-(Vm)2 -f A0( 



(2) 



based on a magnetization like order-parameter m(r). 
A bulk potential 0(m) describes the coexistence of the 
phases a and (3 which, on imposing Ising symmetry, 
we identify with the spontaneous magnetizations —mo 
and -|-mo respectively. The relative potential A0(rn) = 
(j){m) — (j){mQ) is introduced to remove from the total 
free-energy a contribution proportional to the volume. 

From @, we wish to derive the form of an interfacial 
Hamiltonian pertinent to an "m'*" bulk potential 



(pirn) 



(3) 



2 4 

below the bulk critical temperature. Bulk-like fluctua- 
tions are treated in mean-field fashion so the latter con- 
dition implies r > 0. Written in terms of the mean- 
field spontaneous magnetization mg = \/rJu, and inverse 
bulk correlation length n — ^/2r^ the relative potential is 



A<?;>(m) 



8m^ 



mo) 



(4) 



For later purposes, it is convenient to re-express this as 

6m 1 5m? 
"^4 ^ 



A(/)(m) = — 5m^ 



1 



(5) 



mo "t iii^ 

where we have defined 5m = |m| — mo. A much easier 
starting point for analysis is the DP approximation 

/I 



A(/)(°)(m) = — (|m| - mo) 



(6) 



in which one neglects the higher-order cubic and quartic 
terms Q . The superscript indicates that the LGW model 
with a DP potential is the starting point for our pertur- 
bation theory. We define the reference Hamiltonian 



i7(0)r 



dx 



1 



(Vm)2 + A(?!)(")(m) 



(7) 



Finally, we adopt the simplest choice of boundary condi- 
tion corresponding to fixed surface magnetization. Thus, 
if r^, = (x, V'(x)) denotes an arbitrary point on the wall, 
we require 



m(r^) = mi, 



(8) 



for a fixed value of mi . Without loss of generality, we set 
mi > so that the wall preferentially adsorbs a thin layer 
of net positive (/3-like) magnetization. For both the DP 
and the full "m*" models, the mean- field wetting phase 
boundary, as defined for a planar wall, corresponds to 
mi ~ mo. That is, for mi < mo the planar wall-/? inter- 
face is partially wet (with finite contact angle Q) and the 
a 1/3 interface unbinds continuously as mi m^ . The 
wall is completely wet by the f3 phase for mi > mo cor- 
responding to 6* = 0. Note that, with these boundary 
conditions, first-order and tricritical wetting transitions 
do not arise. It is convienient to introduce a dimension- 
less temperature-like scaling field 



t = 



mi — mo 
mo 



(9) 



which measures the deviation from the three-dimensional 
critical wetting phase boundary. 

Following Fisher and Jin 0, Q we define the interfa- 
cial co-ordinate €(x) using a crossing-criterion and con- 
sider constrained profiles which have a surface of iso- 
magnetization m-^ = at some prescribed interfacial 
configuration: 



m{rt) = 



(10) 
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for all points ri = (x,f(x)). Starting from a suit- 
able microscopic model Him], the interfacial Hamilto- 
nian H[£,tp] is defined via a partial trace over this class 
of constrained profiles. A saddle point approximation 
leads to the Fisher-Jin identification 



(11) 



where F^p [tp] is the excess-free energy of the wall-/? inter- 
face. Here ms(r) is the constrained profile that minimises 
the LGW Hamiltonian with the appropriate boundary 
conditions at the interface, wall and in the bulk. That is, 
the constrained profile satisfies the variational equation 



6H\m] 



(12) 



We now focus on the properties of the the DP model ([7]) 
summarising the main results of our previous article. The 
variational equation leads to the Helmholz equation 



(0) 



(13) 



where Sm^ = \ms.\ ~ vtiq. Again, the superscript (not 
used in 1]) serves to indicate that the DP potential is 
the zeroth-order term in a perturbative expansion. The 
linearity of these equations simplifies considerably the 
derivation of the non-local model. Making use of the 
divergence theorem, 

i7(0)[4«)].-^ /ds,, Vmi^^n^ (14) 

mo /■ (0) "^0 /■ , „ (0) 

2 Je- " 2 Ji+ 

contains only surface terms. Here, Snii — mi — mp while 
and denote the (local) unit normals at the wall 
and interface respectively, pointing into the bulk. Simi- 
larly, the infinitesimals ds^ and ds£ represent local area 
elements at the wall and interface. The solutions to the 
Helmholtz equations in the bulk (m^ < 0) and wetting 
(ms > 0) regions are written in terms of the Green's 
function 



27r ri - r2 



-ft |ri-r2| 



(15) 



which satisfies the Ornstein-Zernike like equation 



i-Wl^+K^) K{ri,r2) = 2K5(ri-r2) (16) 

and decays to zero as |ri — — > oo. We represent this 
Green's function diagrammatically by a straight thick 
line with the open circles denoting the end points 



Using the Green's function, we identify the constrained 
magnetization in the bulk region 



(0) 

m~ = —mo 



mo 



(18) 



and, within the wetting layer, via the expansion 



dm-"^ — —mo 



-Srrii 



(19) 



In this diagrammatic notation, the wavy lines represent 
the constrained interfacial configuration (top) and wall 
(bottom) , while a black dot on a surface means one must 
integrate over all points on that surface with the appro- 
priate infinitesmal area element. These expressions are 
exact solutions to the Helmholtz equations, and satisfy 
the boundary conditions at the interface and wall to ex- 
ponentially accurate order in the radii of curvature. 

After substituting into (fT4|) and making use of the 
method of images, we arrive at the desired result 



where 



(20) 



(21) 



is the interfacial Hamiltonian of the free a\(3 interface, 

AajS is the interfacial area and Ylf^ = kitiq is the DP 
result for the interfacial tension. The binding potential 
functional is given by 

oo 

W^'^ [£, ^] = E { + ^^"^ + ^2°^ ""+1 } 

n=l 

(22) 

with geometry independent coefficients 



,(0) 



do) 



2i, 



1, 



(23) 



and correspond to surface integrals over products of 
the kernel K. The diagrammatic and algebraic represen- 
tations of the first three terms are 



ds^dsi K{r.^,re) (24) 



I 



(17) Vl\%^\=^^[^ 




(25) 
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ds^ K{r^,ri) )■ (26) 



although for certain configurations further simphfication 
is possible (see later). All diagrams have a zig-zag struc- 
ture, for example, 



The coefScients appearing in these expressions are deter- 
mined by the coefficients of the binding potential func- 
tional. Thus, we identify w 



w. 



(0) 
20 



and s 



(0) 
21 



(0) 
10 



^(0) 
'10 



,(0) 



a\ ' while 



2fei°''. We now wish to see 



how the above results are altered when one goes beyond 
the DP approximation. 



III. FEYNMAN-HELLMAN THEOREM AND 
PERTURBATION THEORY 



Thus, up to "two tubes" , the asymptotic expansion of W 



,(0) 



•(27) 



with the higher-order terms resumming to give the hard- 
wall repulsion of the wall. 

The coefficient of each diagram ilj^ has the dimensions 
of a surface tension and exhibits a power-law dependence 
on the scaling field t. A handy rule to remember is that 
the power of t is the same as the number of black dots 
on the wall that are singly connected, i.e. have only one 
kernel K attached to it. This rule also applies to all dia- 
grams containing only black dots which will be generated 
in the perturbation series described next. 

One also generates an expression for the excess free- 
energy of the wall-/? interface 



T^ + i (28) 



R 



which depends on the area Ayj and the mean-curvature 
of the wall. The latter are expressed in terms of the local 
principle radii of curvature and r1^ . The tension and 
rigidity are given by — Kva^t^ /2 and C^°^ = 
respectively. 

Finally, for wetting at planar substrates [ip = 0) the 
non-local model recovers the known form of the approx- 
imate local interfacial Hamiltonian when V€ ^ 1. We 
find 



ij(o)r 



ri](o)(£) 



(V^)2 + VFi")(f) 



y(0) A 



(29) 

where W^(^) and S(^) = E^^ + AS(^) are the binding 
potential function and effective position-dependent stiff- 
ness, respectively. Within the DP approximation. 



and 



w 



(0) 
10 



■ w. 



(0) 
20 ' 



(0) - 

Sj^Q e 



JO) 
'21 



-2K.I 



(30) 



(31) 



and are identical to the findings of Fisher and Jin 0, 0] , 
who derived the small gradient (local) limit directly. 



Let us suppose that our microscopic model H[m] can 
be written 



(32) 



containing a dimensionless field e which will later act as 
a small parameter. The reference Hamiltonian is the DP 
model, while 



H^^\m] = jdr A$(i)(m) 



(33) 



accounts for cubic and quartic corrections obtained by 
writing 



A(/)(to) = 



{\m\ — moY 



-HeA$(i)(m) (34) 



with 



A$(i)(m) 



■ Sm 



Sm 
Too 



1 Sm 



4 TOq 



(35) 



Thus, the potential (I34p interpolates between the DP 
model (e = 0) and the "m"'" model (e — 1). Recall that 
the interfacial model is identified by evaluating H[m] at 
the constrained magnetization m^. Taking the derivative 
of the constrained Hamiltonian 



^ H^^>[ms]+ dr — 

de J dm 



dm-= 



(36) 



which, by virtue of the variational condition (I12D . leads 
to the familiar expression 



dH[m^] 



= Idv A^'^^\m-E), 



(37) 



similar to the well known Feynman-Hellman theorem in 
standard quantum mechanics. Note that the functional 
on the R.H.S depends on the full (e dependent) con- 
strained magnetization and is a convenient means of for- 
mulating a perturbation expansion 



^2^(2) 



(38) 



From this, it is straightforward to determine the corre- 
sponding expansion for the binding potential functional 

W[l, V] = V^] + e W^^^ [1, 1/;] + M/(2) [£, i/;] -f ... 

(39) 
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where the leading order-term is the DP result ((22)1 . In 
addition, we will also be able to compute expansions for 
the free interface Hap [I] and the excess free-energy of the 
wall-/? interface 

To determine the first-order and second-order pertur- 
bation functionals i/*-^-* and H^^'' , we return to the Euler- 
Lagrange equation for the constrained profile 



c)A$(i)(mH) 



dm 



and seek a perturbative solution 

6ms.{r;e) = 6'm~\r) + e Sm~\r) + 



(40) 



(41) 



By definition, the leading-order term is the DP result, 
which satisfies the Helmholtz equation (fT3|) . while the 
first-order correction satisfies the inhomogeneous PDE 



(1) 



(1) 



9A$(i)(TOi''^) 



dm 



(42) 



and vanishes at the interface, wall, and at infinity. Com- 
bining these, we obtain 



= Jdr A$(i)(mi''^) 



and 



dr 6m 



1) 9A$(i)(? 



,(0)n 



dr 



(43) 



(44) 



A simplifying feature of the first-order correction is that 
it only depends on the zeroth-order profile m-'^ as calcu- 
lated within the DP approximation. We begin with such 
a calculation for some preliminary quantities. 



IV. FIRST-ORDER PERTURBATION THEORY 
FOR THE FREE HAMILTONIAN 

Consider a free but constrained configuration of the 
a — (3 interface. That is, the interface is infinitely far 
from any confining wall but the magnetization is con- 
strained to be zero along a surface at height ^(x). Bulk 
phases a and (3 lie above and below the interface re- 
spectively. The interface partitions the system into two 
regions whose order-parameter fluctuations are shielded 
from each other, by virtue of the crossing-criterion. The 
zcroth-ordcr DP expressions for the position-dependent 
magnetizations in these regions are 



m~^ (r) = —mo + m^ 



and 



(0)/ \ 
ml (r) 



mo - mo 



(45) 



(46) 



above and below the interface respectively. The first- 
order result for the free interfacial Hamiltonian is 

[^] = dr A$(i) (mP )+ ■■■ (47) 

where the first term is simply the zeroth-order DP result 
= 4°/3^"/5. Hence, 

where we have expressed the results diagrammatically. 
The single wavy line represents the free interface while 
the thick straight lines denote the Green's function K. 
The diagrams appearing in this formula are all of the 
same type and have n = 3,4 (black) dots on the in- 
terface and one (black) dot either above or below the 
surface. They correspond to multi-dimensional integrals. 
For example. 



^kJ dr S^Jdsi K{ri,r)^ 



(49) 



where, in general, the integrand contains n Kernels K 
connecting a point off the interface to n different points 
on it. Black dots on the surface have the same intepreta- 
tion as before - one must integrate over all points on the 
surface with the appropriate area element. A black point 
off the surface means that one must integrate over the 
appropriate semi- volume V+ (here above the interface) to- 
gether with a multiplicative factor of k. The latter is 
introduced so the diagram has the dimensions of area. 
Again, each Kernel may be thought of as representing a 
short tube-like fluctuation protruding from the surface, 
only a few bulk correlation lengths long (since the Kernel 
decays exponentially quickly). Such fluctuations can be 
thought of as giving the interface a " corona" . As we shall 
show, these shift the DP expression for the surface ten- 
sion and also introduce curvature corrections. To see this, 
consider first the case of a planar interface of (infinite) 
area Aap. By definition, the value of the Hamiltonian 
per unit area is equal to the surface tension, so we can 
identify 



The integrals are easily performed 



which implies the tension is shifted to 

1 1 
16 



Sa/3(e) 



1 



(50) 



(51) 



(52) 
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where we have highhghted the different numerical contri- 
butions for the cubic and quartic perturbations. Set- 
ting e — 1, we find « 0.72 Kmg, which is in 
much better agreement with the mean-field expression 
= (2/3)Km^ of the fuU "m^" theory [l^. Thus, the 
dominant numerical correction to the DP expression for 
the surface tension arises from the cubic term in Ai/)^^) 
and is accurately accounted for by first-order pertuba- 
tion theory. This point is further amplified by calculat- 
ing exactly the mean-field surface tension SQ^(e) for the 
potential 



SQ/3(e) 



4(1 



In- 



4 

2(1 + 



2 (V4 - 3e - 2) 



(53) 



described for the free interface and, to first-order, we find 



J +3 



(55) 

where this time the wavy line denotes the shape of the 
bounding wall. The diagrams are easily evaluated as an 
expansion in the inverse principal radii of curvature at 
the wall, and we find 




(56) 



ft is straightforward to check that this is consistent with 
the limiting cases at e = 1 and e = respectively, and 
also reproduces the perturbation expansion (|52p . While 
this function looks rather ominous, it is almost linear in 
character over the required domain. 

In addition to correcting the value of the surface ten- 
sion, the "corona" diagrams lead to curvature terms, 
which reveal the more general structure of the free Hamil- 
tonian. To appreciate this, consider the case of an un- 
dulating interfacial profile. Provided the local principal 
radii of curvature i?f (x), i?|(x), are always large, one can 
expand the integrals to find 

(54) 

where Kap = eTO§/64K, and — — emQ/128K are 

the bending rigidity and saddle-splay coefficients of the 
square mean- curvature and Gaussian curvature, respec- 
tively [1^ |2l[. The notation here is similar to that 
adopted by Blokhuis and Bedeaux [2l[ , although we have 
added a subscript to try to avoid confusion with the in- 
verse bulk correlation length. Note there is no term pro- 
portional to the mean-curvature as required by the Ising 
symmetry. 

A similar calculation reveals the general structure of 
the wall-/3 interfacial free-energy. Consider the interface 
between a wall described by the height function 'ipi'^) and 
the bulk (3 phase corresponding to spontaneous magneti- 
zation Too- Recall that the magnetization at the surface 
TOi is positive so that this interface does not exhibit any 
wetting behaviour. The DP result, Eq. (pS)) . for the ex- 
cess free-energy involves only the area and local mean 
curvature of the wall. No higher order curvature cor- 
rections are present. Beyond DP approximation, we may 
reasonably expect this to change with the cubic and quar- 
tic interactions giving rise to additional curvature contri- 
butions. The perturbation theory is very similar to that 



where the ellipses denote higher-order terms in the curva- 
ture. The new surface tension Y,^p and bending rigidity 
coefficient C^jp contain very small corrections of order 
0{et^) to the DP results quoted earlier. The new rigidi- 
ties K^p ~ Kwji are 0{tt^) and are considerably smaller 
in magnitude than Cwp- 



V. FIRST-ORDER PERTURBATION THEORY 

FOR W[t,ip] 

A. General equations 

To begin, we restate the perturbation theory for the 
bulk potential in a slightly more general way. The calcu- 
lation of the non-local binding potential only requires us 
to specify the form of the bulk potential in the wetting 
layer region where to > 0. We write 

n—3 

where Sm — m — toq and the e„ are all dimensionless 
parameters. Thus, the usual "to*" theory corresponds 
to £3 = ^,£4 = I and e„ = for n > 4. To first-order 
in perturbation theory, all the contributions are additive 
and we seek to write the non-local binding potential func- 
tional 

00 

W[e, tp] = [i, i>] + J2^n Wi^^ [£,tp]+ ■■■ (58) 

n=3 

where, in an obvious notation, the Wn^^ are the pertur- 
bations corresponding to the term Sm" in the bulk po- 
tential. To determine these, it is convienient to order the 
expansion of Sm-"^ in the number of tubes that span the 
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interfaces 



and 



Sm-^ — I Smi 



mo 



(59) 



From it follows that the first-order perturbations 

are given by 

M/^i) [£, ^] = ml'- J dr (<5m^°)) " - [£] - [^j] 

(60) 

where V/j denotes the volume of the wetting layer between 
the wall and interface. The functionals A[i] and B[ip] 
do not describe interactions between the interface and 
wall and are introduced so that W vanishes for infinite 
separation. For example, 



(61) 



(62) 



where, in each case, the wavy line denotes a configuration 
of the surface that corresponds to the argument of the 
functional. All that remains now is the evaluation of the 
integrals and the classification and simplification of the 
ensuing wetting diagrams. 



B. Wetting diagrams for Cubic and Quartic 
interactions 

Substituting the magnetization profile into the first- 
order perturbation expression (|60p for n — 3 and n — A 
leads to the following expressions for the first-order cubic 
and quartic corrections to the DP functional: 



(1) 



(1) 



(64) 



6t' 




-6t' 



Higher-order diagrams exist but involve at least three 
tubes that span the surfaces and would generate terms of 
order 0{e~^'^^) in the standard binding potential func- 
tion. Each of the new wetting diagrams has one black 
dot lying between the surfaces and represents an integral 
over the volume Vp of the wetting layer. The associated 
infinitcsmal measure is ndr. Thus, the first wetting dia- 
gram in the expansion of ^'3^'' is 



K Jds^ J dr K{r^,r) | Jdsi K{r, ] 



(65) 



where we have labelled the points in an obvious notation. 
It is natural to interpret this as a splitting of a tube-like 
fluctuation connecting the surfaces. The second diagram 
in the same cubic interaction does not involve a splitting 
but instead adds a " corona" corresponding to short tube- 
like fluctuations away from the interface: 

^ = (66) 



(63) 



-^3 



3t 



ds^ds'gdr K{r^, r^) K{r'^, r) | Jds^ Kiv, r^) 



Similar intepretations apply to all the wetting diagrams. 
One contribution which is of particular novelty is the X 
diagram 



(67) 



ds^ds'^dsfds[dr K(r^, r) K{rL, r) K{r, r^) K{r, r^) 



and arises from the quartic interaction. This has an ap- 
pealling physical intepretation as a local pinching of two 
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tubes that span the surfaces. As we shah see, this is a 
rather interesting diagram even though ultimately it does 
not influence the leading-order physics. 

C. Wetting diagram relations 

The cubic and quartic interactions appear to give rise 
to a plethora of new wetting diagrams. However, the 
physics represented by these perturbations is rather sim- 
ple and can be elegantly expressed in a more concise fash- 
ion. The essential ingredients in this simplification are 
various relations between the wetting diagrams which ex- 
press their reducibility. We will illustrate this with some 
examples. 

Consider the first wetting diagram appearing in ^3^^. 
To begin, suppose that the wetting layer has planar area 
Au, and is of constant thickness £. The integrals are easily 
evaluated yielding 

Y = A^e-'"' - e-^'^') (68) 

This can be expressed diagrammatically 

X ^ X - ^ 

showing the perturbative diagram is reducible to the DP 
contributions il} and il^. The net effect of this diagram 
is, therefore, to simply shift the coefficients 

aj"' ai — af'' + 1 K uiq (70) 

appearing in the DP expression for W. Moreover, a nice 
feature of the perturbation theory is that there is no need 
to keep precise book-keeping concerning such shifts. This 
can be done exactly at the end of the calculation once the 
general diagrammatic structure has been elucidated. 

The above expression is not quite the whole relation 
for the wetting diagram since interfacial and substrate 
curvature are not allowed for. More generally, one finds 
(after a few integrations) 

X^X4X:-X+ 

where we have introduced a new type of diagram con- 
taining a black triangle. The triangle will always lie on a 
surface and is intepreted as an integral over the surface 
with local measure ds multiplied by the sum of the local 
principal curvatures, measured in units of k (to preserve 
the units of the diagrams). Thus, 



and similarly if a triangle is placed on the wall. The el- 
lipses in the wetting diagram relation ((72|) denote higher 
order curvature terms which are negligible. 

Similarly, for the second wetting diagram in VFg^^ one 
can write the relation 

:k^^x-i^x— <-) 

where here the ellipses also include terms involving four 
tubes that span the surfaces as well as higher-order cur- 
vatures. The same process is also valid for diagrams with 
two tubes spanning the surfaces. For example 



Again the effect of these diagrams is to shift the coef- 
ficient of the fli diagram and add negligible curvature 
terms. In the first-order perturbation theory all bar one 
diagram can be recast as a sum of the DP diagrams fi}, 
ill s-^d ^2 together with curvature corrections. The only 
contribution for which there is no such relation is the X 
diagram describing the two-tube pinching process (j68p 
which is not reducible. However, relations involving it 
do emerge at second-order in perturbation theory. 

In summary, three effects emerge at first-order in per- 
turbation theory: 1) Rescaling of the coefficients ai, bi, 
etc. 2) appearance of curvature corrections and 3) intro- 
duction of non-zig-zag diagrams describing tube interac- 
tions. 



VI. SECOND-ORDER PERTURBATION 
THEORY FOR W 

At second-order, by far the most important contribu- 
tion arises from the cubic interaction in AC''^^^. Con- 
tributions of order e|, as well as mixing terms £364, are 
small and do not introduce any new physics. For ease of 
presentation, we suppose that the potential perturbation 
has only one power, A<I>(^) = ttiq"" Jm", and deter- 
mine the second-order term in 

W[e, 1/;] = I¥(") [i, V] + e« [e, ^] + el 1^(2) ^] (7g) 

Setting n — 3 at the end of the calculation reveals the 
dominance of the cubic interaction at this order. The 
second-order perturbation is 

V;] = Hjf^^ Jdr 5m^i\dm'i^r-' -~Ai;?[e]~Bl^- 

(77) 

where, as in the first-order perturbation theory, function- 

(2,) (2) 
als A'n'ie] and Bk'Oj] are introduced so that, by con- 
struction, W'^^^ vanishes when the interface is delocalised 
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from the wall. They need not be specified explicitly, as 
they are automatically generated by the integral in ([77)1 . 

The second-order term in the potential W depends on 
the first-order correction to the profile Sm~^ which sat- 
isfies 



V'Sr 



.(1) 



2 c- {1) I 2 'Z — n / <• lu J \ 

— K Om~ +nK mf^ (dm- ) 



2 ™2-n , 



A0)\n-1 



(78) 



Substitution of the DP profile leads to the PDE 



1 n-1 1 n-2 



1 n-1 , 1 n-1 1 n-2 > .. 

where we have curtailed the expansion at two tubes span- 
ning the surfaces, and neglected terms of 0(t'^). The in- 
homogeneous PDE can be solved in a standard manner 
using the same Green's function K{ri,r2). Thus, the so- 
lution can also be written diagrammatically and, after 
some algebra, we find 



/ 1 n-1 



In 1 n-1 



In 1 n-1 



1 n-1 1 n-2 



Specialising in the dominant cubic interaction [n — 3), 
we find for the second-order perturbation in W: 



{4:tDl+ D^^} +0{t^ 



(81) 



where the Dl and denote the following sum of dia- 



grams: 



(82) 



and 



(83) 



(84) 



These diagrams determine the rescaling of the coeffi- 
cients Oi and &i, and also generate curvature corrections 
due to the interface. Again, the key to understanding 
their net effect is through wetting diagram relations. For 
example, the following quartic diagram can be expressed 



(85) 



where the ellipses include higher-order interfacial curva- 
ture terms and four-tube diagrams. In this way, each of 
the contributions in (|8ip can be written as a sum of the 
diagrams 



(86) 



similar to the first-order perturbation theory. If one ex- 
tends the calculation to allow terms of order t^, t^, etc. , 
one also encounters wetting diagrams where corona-like 
tubes eminate from the substrate. These are, in fact, the 
same as the diagrams in and Df but with the interfa- 
cial and substrate surfaces switched. Thus, for example. 



(80) 



(87) 



has a coefficient proportional to and will add higher- 
order powers of t in the expansion of ai, and also generate 
curvature corrections due to the substrate which can be 
recast in terms of the diagram 



1 



1 



1 



(88) 



Again, the general structure obtained from the first-order 
perturbation theory is unchanged. 

Working to 0{t^), one also generates wetting diagrams 
which are closely related to the two-tube pinching pro- 
cess which arose in the first-order perturbation from the 
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quartic interaction. For example, 



(89) 



whose coefBcient is proportional to t^. The two cen- 
tral black dots in the wetting layer are connected by a 
tube-like fluctuation which does not attach to either the 
wall or interface. The connecting tube is necessarily of 
short length because the corresponding integral is heav- 
ily damped by the Kernel K. This is neatly expressed 
diagrammatically 



(90) 



leading to the rescaling of the coefficients of Vl\ and 
X . Curvature corrections, represented by the ellipsis, are 
of negligible importance for two-tube diagrams. 

In summary, second-order perturbation theory leads to 
the same three effects highlighted in the first-order calcu- 
lation: the rescaling of coefficients, and the appearance 
of curvature and tube-interaction diagrams. 



VII. THE GENERAL STRUCTURE OF THE 
NON-LOCAL BINDING POTENTIAL. 

The general structure of the non-local binding poten- 
tial functional for short-ranged wetting is now apparent. 
Up to "two tubes" , the functional has an asymptotic large 
distance decay described by the diagrams 



^T^ + Cl -h C2 ^T^ 



(91) 



which should be compared with the DP result (|27|) . Thus, 
going beyond DP generates curvature terms (shown for 
one-tube diagrams only) and also a tube-interaction dia- 
gram. The corresponding algebraic expressions are given 
by equations (I24|,(l73l),(l88]),(l25|,(l26l) and ([M^ The co- 
efficients are geometry independent and all have power 
series expansions in the scaling field t. The leading-order 
behaviours are 



ai 



Cl 

Km? 



at. 



It, 



C2 



/3, 



di 



Xt' (92) 



and are specified by just four dimensionless constants re- 
flecting the surface exchange symmetry of W. The coef- 
ficients b2, C2 and c?i all vanish as implying that the as- 
sociated diagrams are of negligible importance at critical 



wetting. The second diagram, describing the curvature 
correction due to the a\(3 interface, is necessarily much 
smaller than 17} and is therefore also negligible given that 
Cl also vanishes at the critical wetting phase boundary. 
Thus, the diagrammatic expression for W is the same as 
calculated using the DP approximation in [l| but with 
different numerical coefficients. This is the main finding 
of our study. 

The exact values of the above coefficients can be calcu- 
lated for the "m^" LGW potential ([U, by matching with 
mean-field results for specific interfacial and wall config- 
urations. Consider for example the simplest situation of 
a fiat wall, ip = and a flat interface ^(x) = i. The cor- 
responding planar binding potential function is defined 
as 



w[e, 0] 



A,, 



(93) 



and can be identified with the diagrams 



A^WAi) = ai 



(94) 



The first three diagrams are of DP-type and were dis- 
cussed in [l|. The new diagram can also be evaluated 
exactly 



X 



= A, 



(95) 



implying that there are non-purely exponential terms in 
the binding potential. Thus, the binding potential func- 
tion necessarily has the general expansion 



W^{£) = aie 



-Ki 



diK£)e 



-2kI 



(96) 



with coefficients specified in (|92|) . This is identical to 
the findings of Fisher and Jin who calculated W^r {£) di- 
rectly (5,]. One advantage of the Green's function ap- 
proach is that the diagram leading to the non-purely ex- 
ponential term is isolated and can be evaluated for other 
geometries. For example, for spherical interfacial and 
wall shapes 



(97) 



where, as in [H, A^ = 47ri?2 and A^p = Att {R + £f 
are the areas of the wall and interfacial configurations, 
respectively. 

We can now determine the coefficients ai , bi, ... by 
comparing (|96p with the known asymptotic decay of Wt^ 
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for arbitrary potentials A<I>(m). This can be calculated 
independently without recourse to perturbation theory. 
For planar interfacial and wall configurations, the con- 
strained profile ms = m^(z;£) satisfies the "energy- 
conservation" condition 



dz 



and 



/«x/rT(W(^ + ^).-~' (104) 



= J dxyrrrWe"'"^ (105) 



This can be integrated, and the large distance expansion 
exactly determined. For the "m*" potential we find 



ai 



At, 



bi 



di 



6t' 



(99) 



The curvature coefficient 7 = —8 can be determined in a 
similar fashion by considering spherical wall and interfa- 
cial configurations. 

One can go further in this analysis and determine the 
coefficients for the perturbative potential p4)) to all or- 
ders in e. We only quote the results for ai and bi 



ai 



8t 



bi 



16 



2-e + V4^^ 



(2 



(100) 

which smoothly interpolate between the DP and "m^*" 
theory results. 



VIII. NON-LOCAL MODEL FOR WETTING AT 
PLANAR SUBSTRATES 

In this section, we show how all the wetting dia- 
grams appearing in the asymptotic expansion (j9ip sim- 
plify when the substrate is planar (-0 = 0). Clearly, there 
is no contribution from substrate curvature and we write 
the interfacial model 

H[e]^ H^p[e] + w[e] (loi) 

with planar binding potential functional {W[i] = W[£, 0]) 



W[e] = ai 



•(102) 



containing two new diagrams compared to the corre- 
sponding DP expression. Three of these diagrams can 
be evaluated by simply holding the dot (or triangle) on 
the upper interface fixed and integrating over the wall: 



(103) 



which are all local contributions to the effective Hamilto- 
nian H[£]. In particular, if Vi is small, one can see how 
each contribute towards a local binding potential func- 
tion and/or effective position-dependent stiffness. Note 
that, if the coefficient C2 of the curvature diagram (|104p 
is zero, as it is at DP level, the V,l diagram p03p deter- 
mines the leading order exponential decays of Wj^ {£) and 
S(i?). Beyond DP, however, the curvature diagram also 
contributes to S(^). 

In contrast, the two remaining diagrams, fi^ and X, 
are strongly non-local. As remarked in Q application of 
the convolution theorem reduces the triple integral ([25]) 
to a double integral 



= dsids2 e-"^("^)5(xi2;li2) e-"'("^) (106) 



where £12 = (^(xi) + £{x.2))/2 is the mean interfacial 
height of the two points at the interface. Here S' is a 
two-body interfacial interaction which decays as a two- 
dimensional Gaussian 



12 



A£ 



(107) 



and which controls the repulsion of the interface from 
the wall. By construction, the integrated strength of 
S is unity. There are two features about this effective 
many-body interaction which are worth commenting on. 
Firstly, its range increases as the square-root of the film 
thickness and, therefore, becomes longer ranged as the 
interface unbinds. It is this that necessitates a non-local 
treatment of short-ranged critical wetting, and is respon- 
sible for the breakdown of local theories. Also, the same 
Gaussian interaction (|107p follows from a simple saddle- 
point evaluation of the integral (pS)) over the wall. This 
means that the interaction between two fixed points on 
the interface arises due to a connecting tube that refiects 
off the wall and is of minimal length. This physical in- 
terpretation will be useful in discussions of wetting at 
non-planar walls, where an exact convolution evaluation 
of ^\ is not available. 

Similar arguments apply to the X diagram, describing 
the two-tube interaction which can be written 



= jjds^ds2 e-«^(^^)x(xi2;42)e-^ 



-Kf(x2) 



(108) 



12 



The two-body interaction describing this interaction also 
depends on the mean-interfacial height only, and is given 

by 



(109) 



where r(0, z) is the incomplete gamma function. At 
large distances, this decays similar to the two dimen- 
sional Gaussian p07p . 

Finally, we mention that in the strict small gradient 
limit the non-local Hamiltonian reduces to 



m 



{Vif + WAi)} +^c.t3A^ (110) 



where the position dependent contributions to the bind- 
ing potential and stiffness coefhcient have the general de- 
cays 



Wi£) = wio e-'' + {w2iKe + W2o) e"^"' + ■■■ (111) 



and 



AE(^) 



sio e 



S22 K^f + S21 kI + S2o) e"^'"" + 



(112) 

respectively. All seven coefficients exhibit power-law 
dependences on the scaling field t and are determined 
by the five coefficients ai, bi, ci, 62 and di. We find 
wio ^ sio ~ ti W21 ~ S22 ^ and all other coefficients 
finite at t = 0. These are in precise agreement with the 
local theory of Jin and Fisher 



IX. DISCUSSION 

In this paper, we have extended our earlier derivation 
of a non-local interfacial Hamiltonian for short-ranged 
wetting beyond the DP approximation. We have shown 
that the diagrammatic method introduced in our first 
paper combines rather nicely with perturbation theory 
and allows us to derive the general structure of the bind- 
ing potential functional W. While this contains some 
new diagrams describing curvature corrections and tube- 
interations, these are of negligible importance for wet- 
ting. The dominant diagrammatic structure of W is 
the same as that derived for the DP model, albeit with 
slightly shifted coefficients. The values of these coeffi- 
cients have been determined exactly. We also showed how 
all the diagrams in the asymptotic expansion of W (up 
to two tubes) simplify for wetting at planar substrates. 
While some local contributions are accurately described 
by a binding potential function and position dependent 
stiffness, the non-local contributions generate weak long- 
ranged two-body interfacial interactions, which play a 
crucial role at critical wetting. 

To finish our article, we make a number of remarks 
concerning the interpretation, limitations, and further 
extension of the present approach: 



1. - New features. Going beyond DP approximation 
alters the binding potential functional in three different 
ways: rescaling of the coefficients oi, bi . . . , and the in- 
troduction of curvature and tube interaction diagrams. 
These effects could have been anticipated on very gen- 
eral grounds. Even for a free Hamiltonian, going beyond 
DP alters the surface tension and introduces Helfrich- 
like rigidity terms, consistent with morphological expec- 
tations. The first two effects merely mirror those in W. 
Also, there must be a mechanism which generates non 
purely exponential terms in the binding potential func- 
tion W^(^) beyond DP. This is fulfilled by the two-tube 
interaction diagram X. 

2. - Generalisation of the LGW model. One could con- 
sider a slightly generalised LGW model in which the co- 
efficient of (Vm)^ is a function of the order-parameter 
[m . By expanding this function about toq , it is straight- 
forward to show that the corresponding first-order per- 
turbation corrections to W are equivalent to cubic, quar- 
tic, etc. corrections to the DP potential, which we have 
considered explicitly. The same is also true of square- 
Laplacian contributions to H[m] which are equivalent to 
a shift in the value of n plus the aformentioned cubic, 
quartic and higher-order corrections. In both cases, the 
diagrammatic structure of W remains unaltered to first- 
order. Of course, this robustness is to be expected. As 
remarked at the end of our first paper, the form of W is 
necessitated by exact sum-rule requirements. 

3. - Coupling to a surface field and enhancement. While 
we have restricted ourselves to fixed boundary conditions 
at the wall, it is straightforward to extend the Green's 
function method to different kinds of surface boundary 
conditions. In particular, one may consider LGW models 
of the form [l3| 



H[m] = Jdr I i(Vm)2 -|- A0(m) | + Jds^ 0i(m(r^)) 

(113) 

where (/'i(m) = — — him describes the coupling to 
a surface field hi and enhancement g. In this case, the 
divergence theorem generalises eq. (I14p to 

ij(o)[mi°^] f ds^ (mi°^ - mo) 0'i(mi°^) (114) 

2 Ji- 2 Ji+ 

where we have used the appropriate boundary condition 
in the integral over the wall. Similarly, the diagrammatic 
expansion for the constrained magnetisation reads, in DP 
approximation, 



Sm. 



(0) 



-mo 



+ 



(115) 
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where the black squares on the wall denotes a convolu- 
tion of the Green's function K with an auxiliary func- 
tion /i(r^). This accounts for the variation of the surface 
magnetisation, and is introduced to satisfy the boundary 
condition at the wall. For example, for a planar wall, /i 
satisfies the integral equation 



K /i(r^) + 2k mo 



(116) 



which can be solved via a Fourier transform. From here, 
the calculation proceeds as described in jj, within DP 
approximation, and similar to that described herein using 
perturbation theory. The final result for the two domi- 
nant diagrams in W remains unchanged from eq. ([1]) but 
with coefficients 



bi 



K 8 (/iiTOo + giTT-o) 



9 + 



16 ktoq 



(117) 



n-g (2-e + V4~3i)2 



where we have expressed the results as appropriate for 
the potential (|34p . As anticipated, the diagrammatic 
form is very similar to that described for fixed boundary 
conditions but with the advantage that one can now dis- 
cuss first-order wetting {g > —k) and tricritical wetting 
(g = —k). The problem of tricritical wetting is particu- 
larly interesting because the coefficient &i of the second 
diagram, 17^, in ([1]) vanishes. The repulsion from the wall 
is then determined by a diagram which generates a term 
of order e ~^ in the usual binding potential function. 
This corresponds to the next diagram in the series ((9T|) 
and involves three tubes, as is discussed next. 

4.- Tie dominant three-tube diagram. There are sev- 
eral diagrams involving three tubes which contribute to- 
wards the coefficient of e~^'^^ in Wt^II). However, there 
is only one such diagram whose coefficient does not van- 
ish at t = and is, therefore, necessary for the discussion 
of tricriticality. The diagram in question is 



=kJ^ dr I^JdseK{re,r) 



(118) 



and is generated by the cubic interaction in A0(m). 
Here, Vw denotes the volume of the wall. This diagram 
is the next term in the asymptotic expansion (PT|) and is 
strongly non-local. However, analogous to our previous 
discussion of Sl^, the diagram simplifies and generates 
an effective many-body interfacial interaction between 
points at the interface. For example, for a planar wall. 



the integral reduces to 



(119) 



dSi dS2 dS3 e-''"-^'^ e-^"^^'^ e-^"-^'^ r(xi2, X23,Xi3) 



where the three-body interaction 

T{xi2, X23,xi3) cx Ai23 exp(-Ai23 (a;?24 + 2^23^1 + ^13^2)) 



Here, 



Ai 



23 



2(^1^2 +^2i?3+ 44) 



(120) 



(121) 



and we have abbreviated £1 = ^(xi), etc. 

The fluctuation theory of three-dimensional short- 
ranged tricritical wetting requires a renormalisation- 
group treatment of the flow of this three-body interaction 
simflar to that described in 0] for the two-body term 5", 
eq. (|107p . pertinent to critical wetting. We note that 
the range of this three-body interaction also increases as 
the interface depins implying that non-local effects are 
important at this transition. 

5.- Resummation of diagrams within DP. For wetting 
at a planar wall, it is in fact possible to resum all the 
diagrams appearing within the DP approximation, eq. 
([22|) . This, again, makes use of the idea of effective many- 
body interactions, and is possible because in a general 
diagram one can integrate exactly over any flnal black 
dot of the zig-zag that is on the wall. Thus, the DP non- 
local Hamiltonian for wetting in a planar wall, with flxed 
boundary conditions, is given by 



H 



DPI 



where 



dsi {l + te-'^^ f + (122) 

JJdSidS2 /9(xi) Ai"(xi,X2) p(x2) 



p(x) =e-«^^"^ + te 



(123) 



and 

M(Xi,X2) = 5(xi,X2;42) 



(124) 



Jds3 M(xi,X3) 5(X3,X2;42) e-2-^("^) 



and S is deflned in (|107p . The flrst integral generates 
three local contributions: the surface area, and il^. 
At leading order, the total two-body interaction M is 
given by 5, in which case the t independent terms in the 
second integral reduce to fif . In fact, the two-body term 
M is very well approximated by 



M(xi,X2) 



S'(xi,X2;£i2) 



1 



-2k£ 



(125) 
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which shows how the higher-order diagrams resum to give 
a hard-wall repulsion in the two-body interaction. 

6.- Full diagrammatic structure beyond DP. We have 
not attempted to classify diagrams that contain three or 
more tubes. This is much more cumbersome to do be- 
yond the DP approximation and, with the exception of 
tricritical wetting studies, is largely of academic inter- 
est only. Nevertheless, such structure must exist as can 
be seen from the following argument. We have shown 
that the X diagram generates a term K^e"^**^ in the pla- 
nar binding potential W^(^) whose coefficient is oc t^. 
Naively, this suggests that the mean-field excess free- 
energy, obtained by minimising Wtt (i) , contains a higher- 
order logarithmic singularity t^ln|i|. However, such a 
contribution does not exist as can be seen from the full 
mean-field calculation for the LGW model. This must 
mcian that, when evaluated at the equilibrium mean- field 
wetting layer thickness, the term t'^X in exactly can- 



cels with higher-order diagrams which generate terms of 
order tnle^^'^''' and tnie^'^'^^. Indeed, such diagrams can 
be readily identified in the perturbation theory. This is 
strongly suggestive that at least some higher-order dia- 
grams can be grouped together systematically. 

Further extensions and applications of this work, in- 
cluding a discussion of correlation function structure and 
the presence of long-ranged forces, will be discussed in 
future papers. 
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